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^ Abstract 

a: 

Q_i| We systematically investigate constraints on the parameters of the supersymmetric 



type-I seesaw mechanism from the requirement of successful thermal leptogenesis 
in the presence of upper bounds on the reheat temperature Trh of the early Uni- 
verse. To this end, we solve the flavour-dependent Boltzmann equations in the 
MSSM, extended to include reheating. With conservative bounds on Trh, leading 
to mildly constrained scenarios for thermal leptogenesis, compatibility with obser- 
vation can be obtained for extensive new regions of the parameter space, due to 
flavour-dependent effects. On the other hand, focusing on (normal) hierarchical 
light and heavy neutrinos, the hypothesis that there is no CP violation associated 
with the right-handed neutrino sector, and that leptogenesis exclusively arises from 
the CP-violating phases of the C/mns matrix, is only marginally consistent. Taking 
into account stricter bounds on Trh further suggests that (additional) sources of 
CP violation must arise from the right-handed neutrino sector, further implying 
stronger constraints for the right-handed neutrino parameters. 



1 Introduction 



One of the most appealing mechanisms to generate the observed baryon asymmetry of 
the Universe (BAU), ris/n^ ~ (6.10 ± 0.21) x 10~ 10 [lj, is that of baryogenesis via 
leptogenesis. Thermal leptogenesis is an attractive and minimal mechanism, in which a 
lepton asymmetry is dynamically generated, and then converted into a baryon asymmetry 
due to (B + L)-violating sphaleron interactions [2J. The latter exist both in the Standard 
Model (SM) and in its minimal supersymmetric (SUSY) extension, the MSSM. 

The seesaw mechanism [21 H] provides an elegant explanation for the observed small- 
ness of the neutrino masses and, in addition, it offers the possibility of leptogenesis [5]. 
In this case, the lepton asymmetry is generated by the out-of-equilibrium decays of the 
same heavy right-handed neutrinos which are responsible for the suppression of the light 
neutrino masses. In spite of being one of the most simple frameworks where thermal 
leptogenesis can be realised, the seesaw mechanism introduces a large number of new 
parameters (in both its SM and MSSM versions). Although an important amount of 
data has already been collected, many of the seesaw parameters, namely those associated 
with the right-handed neutrino sector, are experimentally unreachable. As discussed by 
many authors, strong constraints on the seesaw parameter space can be imposed from the 
requirement of successful thermal leptogenesis Typically, all these studies relied on 
the so-called flavour-independent (or one-flavour) approximation to thermal leptogenesis. 
In the latter approximation, the baryon asymmetry is calculated by solving a Boltzmann 
equation for the abundance of the lightest right-handed neutrino, and for the total lepton 
asymmetry. Additionally, in the flavour-independent approximation, the only relevant 
CP violation sources are those associated with the right-handed neutrino sector (more 
concretely, the complex i?-matrix angles, when working in the so-called Casas-Ibarra 
parameterisation [7]). 

In recent years, the impact of flavour in thermal leptogenesis has merited increasing 
attention [8]-[l6]. In fact, the one-flavour approximation is only rigorously correct when 
the interactions mediated by the charged lepton Yukawa couplings are out of equilibrium. 
Below a given temperature (e.g. C(10 12 GeV) in the SM), the tau Yukawa coupling comes 
into equilibrium (later followed by the couplings of the muon and electron). Flavour 
effects are then physical and become manifest, not only at the level of the generated 
CP asymmetries, but also regarding the washout processes that destroy the asymmetries 
created for each flavour. In the full computation, the asymmetries in each distinguishable 
flavour are differently washed out, and appear with distinct weights in the final baryon 
asymmetry. 

Flavour-dependent leptogenesis has recently been addressed in detail by several au- 
thors [10]- [HI- In particular, flavour-dependent effects in leptogenesis have been studied, 
and shown to be relevant, in the two right-handed neutrino models [12] as well as in 
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classes of neutrino mass models with three right-handed neutrinos [14J. The quantum 
oscillations/correlations of the asymmetries in lepton flavour space have been included in 
Ref. [10] and the treatment has also been generalised to the MSSM [14J. 

One interesting implication of the flavour-dependent treatment is that in addition to 
the right-handed sector CP violating phases (that is, a complex i?-matrix), low-energy CP 
violating sources, associated with the light neutrino sector, also play a relevant role. Even 
in the absence of CP violation from the right-handed neutrino sector (which would lead 
to a zero baryon asymmetry in the one-flavour approximation), a non- vanishing baryon 
asymmetry can in principle be generated from the CP phases in the Maki-Nakagawa- 
Sakata matrix, ?7mns- Strong connections between the low-energy CP phases of the 
Umns matrix and CP violation for flavour-dependent leptogenesis can either emerge in 
classes of neutrino mass models [H] or under the hypothesis of no CP violation sources 
associated with the right-handed neutrino sector (real R) [131 UB]. In addition, in 
the latter limit, bounds on the masses of the light and heavy neutrinos and on the 
flavour-dependent decay asymmetries have been derived [IE]. The correlation of the 
baryon asymmetry with the effective Majorana mass in neutrinoless double beta (Ov/3/3) 
decays has also been addressed [15]. Another appealing aspect of the flavour-dependent 
treatment is that at least one of the CP sources, the Dirac CP-violating phase, is likely 
to be experimentally measured (one of the Majorana phases could also in principle be 
measured - even though this represents a considerable challenge [T7], while the right- 
handed phases are experimentally unaccessible). 

In the supersymmetric implementation of the seesaw mechanism, further constraints 
on the seesaw parameter space can arise. These are particularly relevant in models of 
local SUSY (i.e. supergravity). First, one should consider cosmological bounds on the 
reheat temperature (Trh) after inflation, associated with the thermal production of grav- 
itinos. Trh has generally a strong impact on thermal leptogenesis, since the production 
of right-handed neutrinos is suppressed if their mass largely exceeds Trh- In fact, due 
these bounds on the reheat temperature, viable thermal leptogenesis will impose strong 
constraints on the seesaw parameter space of locally supersymmetric models. Secondly, 
additional bounds on the SUSY seesaw parameter space arise from low-energy observ- 
ables, namely lepton flavour violating (LFV) muon and tau decays such as lj — > U 7 and 
lj — ► 3 k (i < j) , and charged lepton electric dipole moments. In studies of LFV, thermal 
leptogenesis in the flavour-independent approximation has been discussed in Ref. [T8] . 
and reheating effects have been explicitly included in Ref. [19] . With a potential future 
observation of the sparticle spectrum, the combined constraints on the SUSY seesaw pa- 
rameters from leptogenesis and LFV could lead to interesting information on the heavy 
neutrino sector, which is otherwise unobservable at accelerators. 

In this study, our aim is to investigate the constraints on the parameters of the type-I 
SUSY seesaw mechanism from the requirement of a successful flavour- dependent ther- 
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mal leptogenesis in the presence of upper bounds on the reheat temperature of the early 
Universe. Previous studies of flavour-dependent thermal leptogenesis were conducted in 
the SM, and for a mass range of the lightest right-handed neutrino where only the tau- 
flavour is in thermal equilibrium. In the present analysis, we will work in the context 
of the MSSM extended by three right-handed neutrino superfields. Moreover, for the 
temperatures we will consider, both tau- and muon-flavours are in thermal equilibrium, 
so that in fact all leptonic flavours must be treated separately. We update the Boltzmann 
equations of Ref. [20] to include flavour effects, and point out which regions of the seesaw 
parameter space generically enable optimal efficiency and/or optimal decay asymmetries 
for leptogenesis (focusing on the case of hierarchical light and heavy neutrino masses). 
We then discuss the differences between the flavour-independent approximation and the 
correct flavour-dependent treatment. We encounter interesting new regions of the seesaw 
parameter space, which are now viable due to flavour-dependent effects. On the other 
hand, and as we will discuss throughout this work, scenarios of leptogenesis solely aris- 
ing from the Umns phases are quite difficult to accommodate, and become increasingly 
disfavoured when stronger bounds on Trh are taken into account. 

In the presence of strict bounds on Trh, whether or not it is possible to generate the 
observed baryon asymmetry exclusively from low-energy Dirac and/or Majorana phases 
is still a question that deserves careful consideration. Likewise, it is worth considering 
to which extent right-handed sector phases (other than those of the Umns) could affect a 
potentially viable scenario of low-energy CP violating leptogenesis, and vice-versa. 

Our work is organised as follows. In Section [2} we briefly summarise the most relevant 
aspects of the type-I SUSY seesaw mechanism. Section |3] is devoted to the discussion of 
thermal leptogenesis with reheating. In addition to estimating the baryon asymmetry, we 
focus on the constraints on the reheat temperature arising from the gravitino problem. A 
short summary of the limitations and approximations of our computation is also included. 
In Section HJ we finally discuss the constraints on the seesaw parameters obtained from 
the requirement of successful leptogenesis, namely constraints on the neutrino masses 
and on the i?-matrix mixing angles. Our conclusions are presented in Section [51 

2 The seesaw mechanism and its parameters 

In what follows, we briefly introduce the most relevant features of neutrino mass gener- 
ation via the seesaw mechanism. In the MSSM extended by three right-handed neutrino 
superfields, the relevant terms in the superpotential to describe a type-I SUSY seesaw 
are 
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In the above, N c denotes the additional superfields containing the right-handed neutrinos 
N{ and sneutrinos iVj. The lepton Yukawa couplings and the Majorana mass ttlm 
are 3x3 matrices in lepton flavour space. From now on, we will assume that we are in 
a basis where and itlm are diagonal. After electroweak (EW) symmetry breaking, the 
charged lepton and Dirac neutrino mass matrices can be explicitly written as mi = Xe v\, 
m D = X u v 2 , where Vi are the vacuum expectation values (VEVs) of the neutral Higgs 
scalars, with v%(2) = v cos(sin)/5 and v = 174 GeV. 

The 6x6 neutrino mass matrix, whose eigenvalues are the masses of the six Majorana 
neutrinos, is given by 

y A v v 2 m M J 

In the seesaw limit, ^> v 2 , one obtains three light and three heavy states, Ui and 

Ni, respectively. Block-diagonalisation of the neutrino mass matrix of Eq. ([2]), leads (at 
lowest order in the (X u v 2 /rriM) n expansion) to the standard seesaw equation for the light 
neutrino mass matrix, 

m u = —v\ \ T V rri^ \ v , (3) 

and to m^v — tum- Since we are working in a basis where mu is diagonal, the heavy 
eigenstates are then given by 

m N = m^ ag = diag(m Nl ,m N2 ,m N3 ) . (4) 

The matrix m v can be diagonalised by the unitary matrix ?7mns, leading to the following 
masses for the light physical states 

^MNS^^MNS = «^ mg = diag (771^,771^,771^) . (5) 

Here we will use the standard parameterisation for the £/mns, given by 

Cl2 C13 Sl2 C13 S13 

#MNS = | -«ia C23 - C12 S 23 8 13 e lS C12 C 23 - S 12 S 23 S 13 e lS S 23 C13 j . P , (6) 
s 12$2 3 — C12 C23 Sx3 e tS —Ci 2 S 23 — Si2C2 3 Si3e* 5 C23 C13 

with P = diag (e~ 1 ^, e~ 1 ^, 1), and where Cy = cos%, Sy = sin%. The parameters 0^- 
are the neutrino flavour mixing angles, 5 is the Dirac phase and (px }2 are the Majorana 
phases. 
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In view of the above, the seesaw equation, Eq. (j3J), can be solved for the neutrino 
Yukawa coupling \ v using the Casas-Ibarra parameterisation [7j as 



where R is a generic complex orthogonal 3x3 matrix that encodes the possible extra 
neutrino mixings (associated with the right-handed sector) in addition to those in the 
^mns- R can be parameterised in terms of three complex angles, 9i (i = 1, 2, 3) as 



with Q = cos^j, Si = sin 6*j. Eq. (jTJ) is a convenient means of parameterising our ignorance 
of the full neutrino Yukawa couplings, while at the same time allowing to accommodate 
the experimental data. Notice that it is only valid at the right-handed neutrino scales 
rriM, so that the quantities appearing in Eq. (I7j) are the renormalised ones, m^ iag (uim) 
and LV" MNS (m M ). 

In this study, we shall mainly focus on the simplest scenario, where both heavy and 
light neutrinos are hierarchical, <C mjv 2 <C rn^ 3 and m Ul -C m U2 <C m V3 , and in 
particular, we will assume a normal ordering of the light neutrinos. Thus, the masses m U2 3 
can be written in terms of the lightest mass m Ul and the solar /atmospheric mass-squared 
differences as m 2 U2 = Arn^, + and m^ 3 = Am^ tm + m^. 

In summary, when working in the i?-matrix parameterisation, the 18 parameters of 
the seesaw mechanism are accounted by the three heavy neutrinos masses, mjy-j, the mass 
of the lightest neutrino m Vl plus the two mass squared differences Am^ and Am^ tm , 
the three mixing angles dij and three CP violating phases 8,cpi,cp2 of the Umns matrix, 
and the three complex angles 9i of the matrix R. As mentioned in the Introduction, 
many of the latter parameters, namely those associated with the heavy neutrino sector, 
are experimentally unreachable. Nevertheless, it is possible to derive interesting bounds 
from the requirement of successful thermal leptogenesis, and we proceed to do so in the 
following sections. 

3 Flavour-dependent thermal leptogenesis with 
reheating 

As recently pointed out [TOl [TT1 [12], flavour can have a strong impact in baryogenesis 
via thermal leptogenesis. The effects are manifest not only in the flavour- dependent CP 
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asymmetries, but also in the flavour-dependence of scattering processes in the thermal 
bath, which can destroy a previously produced asymmetry In fact, depending on the 
temperatures at which thermal leptogenesis takes place, and thus on which interactions 
mediated by the charged lepton Yukawa couplings are in thermal equilibrium, flavour- 
dependent effects can have a strong impact on the estimation of the produced baryon 
asymmetry |8J-[16J. For example, in the MSSM, for temperatures between circa (1 + 
tan 2 /3) x 10 5 GeV and (1 + tan 2 /3) x 10 9 GeV, the /i and r Yukawa couplings are in thermal 
equilibrium and all flavours in the Boltzmann equations are to be treated separately. For 
instance, for tan/3 = 30, this applies for temperatures below about 10 12 GeV (and above 
10 8 GeV), a temperature range we will be subsequently considering. Moreover, in the full 
flavour- dependent treatment, lepton asymmetries are generated in each individual lepton 
flavour. Processes which can wash out these asymmetries are also flavour- dependent, i.e. 
the inverse decays from electrons can only destroy the lepton asymmetry in the electron 
flavour. We will address the latter issues in Section 13.11 

In thermal leptogenesis, the population of right-handed neutrinos N\ is produced from 
scattering processes in the thermal bath. To generate the observed baryon asymmetry 
comparatively high temperatures of the early Universe are required, and these should not 
lie much below the mass of the lightest right-handed neutrino, m^. Even under optimal 
conditions, thermal leptogenesis demands > 10 9 GeV (for hierarchical light and 
heavy neutrinos) [21] . High temperatures compatible with thermal leptogenesis can arise 
in the process of reheating which takes place after cosmic inflation [22] . The temperature 
of the Universe at the end of reheating is referred to as the reheat temperature Trh- 
However, particularly in locally supersymmetric theories, T RH is often constrained, as we 
will discuss in Section I3~2l 

In the presence of such bounds on Trh, the requirement of successful thermal lepto- 
genesis imposes severe constraints on the seesaw parameters. In our numerical calcula- 
tions, and following Ref. [20], we include reheating in the flavour-dependent Boltzmann 
equations [HI El EED, EE], generalised to the MSSM [H], in a simplified but compara- 
tively model-independent way. In particular, we assume that the lightest right-handed 
(s)neutrinos are only produced by thermal scatterings during and after reheating. More- 
over, we neglect model-dependent issues such as the production of N\ (and iVj.) during 
preheating, or from the decays of the scalar field responsible for reheating. For complete- 
ness, the technical aspects of the Boltzmann equations are given in Appendix |A] Further 
details on the estimation of the produced baryon asymmetry using the flavour- dependent 
Boltzmann equations can be found in Refs. [HI EH [HI 121 E] • Finally, and concerning 
the inclusion of reheating, we refer the reader to Ref. [2TTj . 

Let us now begin by reviewing (omitting technical aspects) the procedure for esti- 
mating the baryon asymmetry produced by thermal leptogenesis in the MSSM when 
reheating effects are included. 
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3.1 Estimation of the produced baryon asymmetry 



The out-of-equilibrium decays of the heavy right-handed (s)neutrinos give rise to flavour- 
dependent asymmetries in the (s)lepton sector, which are then partly transformed via 
sphaleron conversion into a baryon asymmetry Yg0. The final baryon asymmetry can be 
calculated as [1] 



a 

where Y& a = Y B /3—Y La are the total (particle and sparticle) B/3 — L a asymmetries, with 
Yi a the lepton number densities in the flavour a = e,/x, r. The asymmetries Y& a , which 
are conserved by sphalerons and by the other MSSM interactions, are then calculated 
by solving a set of coupled Boltzmann equations, describing the evolution of the number 
densities as a function of temperature. We consider the simplest case of thermal leptoge- 
nesis, where only the lightest right-handed neutrinos are produced in the thermal batlo 
In the MSSM, the asymmetries for the decay of the lightest right-handed (s)neutrinos 
into (s)leptons of flavour a (defined in Eq. (IA.1O0 ) satisfy Si,a = e i a = e i,a = e ia- Thus, it 
is convenient to write the solutions of the Boltzmann equations in terms of the flavour- 
dependent decay asymmetries £\ >a and flavour-dependent efficiency factors 7] a as 

Y Aa = v a ei,« [^(T>m JVl )+^(T>m Ari )] . (10) 

In the above, Y^ and Y^ q are the number densities of the lightest right-handed neu- 
trino and sneutrino in the Boltzmann approximation (i.e. assuming common phase space 
densities for both fermions and scalars) if they were in thermal equilibrium at T ^> m^, 



Y£(T » m Nl ) » Yg(T » m Nl ) « , (11) 



with = 228.75 denoting the effective number of degrees of freedom. While the equi- 
librium number densities mainly serve as a normalisation, the relevant quantities are the 
decay asymmetries and the efficiency factor, which we now proceed to specify. 



3.1.1 The flavour-dependent decay asymmetries £ lj(X 

In the basis where both the charged lepton and right-handed neutrino mass matrices are 
diagonal, the asymmetries Si >a for the decay of the lightest right-handed neutrino into 

1 Throughout this study Y will always be used for quantities which are normalised to the entropy 
density. 

2 The limitations of this (and other) approximation(s) will be discussed in Section [3. 1.31 
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lepton and Higgs doublets (c.f. Eq. (jA.lOp ) are given by [23J 



ei, c 



1 Sj=2 3 Im [(^f)la (A„ Aj,)u (Aj,) a jl / 

5 



}7T 



(At A,) 
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ni 



A'l 



:i2i 



where 



1 — x 



-In 



1 + x 



X 



2>1 



(13) 



For -C rriN 2 ,mN 3 , and using the i?-parameterisation (see Eq. (JTj)), the decay asym- 
metries in the MSSM can be written as 



3 771^ 



Im 



^2f3p m lp 2 (^MNs)a/j(^MNs)ap -Rl/3 -Rip 

E5 m *>s \ R IS\ 2 



■TTVn 



(14) 



Notice that there is no sum over a in Eq. (TT^|) . which implies that both the £/mns phases 
and a complex i?-matrix can contribute to the CP violation necessary for leptogene- 
sis. This has been recently pointed out by several authors [121 EH E2, [16], and is in 
direct contrast with the flavour-independent approximation, where (working in the R- 
parameterisation) the Umns plays no role in the decay asymmetry E\ = 



3.1.2 The efficiency factors r) a 

As already mentioned, the lepton asymmetries in each individual flavour are obtained by 
solving the set of flavour-dependent Boltzmann equations, taking into account reheating 
effects (c.f. Appendix [A]) . Parameterising the solution of the Boltzmann equations as in 
Eq. ( 1101) implicitly defines the efficiency factor rj. In our approximation, rj is a function 
of the ratio / T^n , of the product A aa m l Q (no sum over a), and of the total washout 
parameter rh\. The quantities mi jCf and rh\ are defined as 

mi, a = (K)ia(*l)ai , rn x = }m lj(X . (15) 

m Nl 

a 

If leptogenesis takes place at temperatures between about (1 + tan 2 /?) x 10 5 GeV and 
(1 + tan 2 (3) x 10 9 GeV, which is the case we will consider in this study, A is given as in 
Ref. |14] (see Appendix [A]). Here we will neglect the small off-diagonal elements of A, 
and use only the leading diagonal entries, 

A » diag(-93/110, -19/30, -19/30) . (16) 
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In addition, there is also a slight dependence of r] on tan (3. We also notice that the 
parameters K and K a (see Eq. (IA.14I) ) are often used instead of rh\ and rh\ i(X . Using the 
i?-parameterisation, the washout parameters rh\^ a can be written as 



The numerical results for the efficiency factor as a function of |^4 aQ |^i, Q and t^a^/Trh are 
shown in Fig. [TJ As can be seen, the efficiency is indeed optimal for values of |^4 aQ |?Tii iQ ~ 
m* , with m* = sin P x 1.58 x 10 -3 eV p3] (see Eqs. (jXUltOll ). and quickly drops for 
either smaller or larger lA^mi^. From Fig. [T], it is apparent that for |^4 aa |mi ja < m*, 
the efficiency exhibits a less pronounced drop if mi/(\A aa \rhi >a ) ^> 1. Moreover, one 
verifies that the efficiency is strongly reduced if 3> Trh; for instance by more than 
three orders of magnitude if w,n x ~ 10 x Trh- This stems from our assumption that Ni 
(and Ni) are exclusively produced from thermal scatterings during and after reheating. 
With respect to Fig. [TJ let us finally notice that even though the curves were obtained 
for an example of tan f3 = 30, the results are approximately the same for other moderate 
(and even large) values of tan/3. Likewise, the contour lines for larger mi/ (\A aa \mi tCt ) 
look virtually the same as for rhi/(\A aa \rhi^ a ) = 100. 

One important difference between the flavour- dependent treatment and the flavour- 
independent approximation is that in the former case the total baryon asymmetry is 
the sum of the distinct individual lepton asymmetries, weighted by the corresponding 
efficiency factor, as in Eq. fllOp . Therefore, the total baryon number is in general not 
proportional to the sum over the individual CP asymmetries £i = r\ ^ 1 CX. 5 *3jS assumed in 
the flavour-independent approximation, where flavour is not considered in the Boltzmann 
equations. Moreover, the flavour-dependent efficiency factors rj a are in general not equal 
to each other, and can strongly differ from the common efficiency factor rj (itself a function 
of the common washout parameter rh\ of the flavour- independent approximation). Taking 
all the latter into account can lead to dramatic differences regarding the estimates of the 
produced baryon asymmetry in models of neutrino masses an mixings p2} [Hj. 

3.1.3 Limitations and approximations of the computation 

As this point, it is pertinent to summarise the approximations and assumptions leading 
to the calculation of the produced baryon asymmetry. 

For the right-handed neutrino sector, we assume a hierarchical mass spectrum, which 
means that we do not take into account the possibility of resonant effects in leptogene- 
sis [21]. Furthermore, we assume that only decays of N\ (and N\) significantly contribute 
to the final asymmetry, thus neglecting contributions from the decays of A^2, (and 
N2,Ns). This is justified when either the washout processes mediated by the lighter 
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logjoflAoJfti^/eV) 

Figure 1: Contour lines for the flavour-dependent efficiency factor r\ in the MSSM as a func- 
tion of |A qq |toi iQ (no sum over a) and toa^/Trh, for fixed values of fhi/(\A ao ,\mi t01 ). We take 
fhi/(\A aa \mi t a) = 2,5 and 100, respectively corresponding to the solid, dotted and dashed contour 
lines. 

right-handed (s)neutrinos are sufficiently active for each flavour, and thus indeed de- 
stroy the asymmetries arising from the decays of the heavier right-handed (s) neutrinos, 
or when T R h <C m n 2 , m n 3 . It is important to notice that under certain conditions the 
heavier right-handed neutrinos can also play a role for leptogenesis [25], and should in 
principle be included. 

We include reheating in the flavour-dependent MSSM Boltzmann equations in a sim- 
plified, but comparatively model-independent way, following Ref. [20]. We notice that 
in specific models of reheating after inflation, the prospects for leptogenesis can be dif- 
ferently]. In particular, we assume that the lightest right-handed (s) neutrinos are only 
produced by thermal scatterings during and after reheating. We also neglect the pos- 
sibility of producing N\ (and Ni) from the decays of the scalar field <fi responsible for 
reheating, or of producing N\ (and N\) during preheating. It is further assumed that 
the maximally reached temperature during the reheating process, T max , as well as the 
mass of 0, are much larger than m^ l . Furthermore, we neglect the potential implications 
of supersymmetric flat directions for reheating and thermal leptogenesis, which are still 
under controversial discussion [28J. 

When solving the Boltzmann equations for flavour-dependent leptogenesis, we focus 
our attention on the case where, during the relevant temperatures for leptogenesis, the 
interactions mediated by each of the charged lepton Yukawa couplings are either fully in 
equilibrium, or out of equilibrium. In the MSSM, for values of around (1 + tan 2 0) x 

3 Constraints on the seesaw parameters differ significantly if, for instance, the inflaton predominantly 
decays into right-handed (s)neutrinos [26j . or if one works in scenarios where the inflaton is one of the 
right-handed (singlet) sneutrinos [57] , 
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10 9 GeV, the reactions induced by the muon Yukawa coupling are close to equilibrium and 
the quantum oscillations of the asymmetries may not have been dumped fast enough to be 
neglected. To take these effects into account, the Boltzmann equations may be generalised 
so to include quantum oscillations [TD]. In this study, we chose tan/? sufficiently large 
so that we can safely assume that the charged /i and r Yukawa couplings are in thermal 
equilibrium, and that all flavours in the Boltzmann equations can be treated separately. 
Furthermore, we neglect the small off-diagonal elements of the matrix A, which appears 
in the washout terms of the Boltzmann equations for Y& a (c.f. Appendix [A]) . We note 
that this approximation is crucial if one wants to introduce an efficiency factor 77, which 
is a function of the ratio m^/T^H, of the product A aa mi jQ (no sum over a), and of the 
total washout parameter rh\. We have numerically verified that in the regions of interest 
for this study, including the off-diagonal elements has only a small effect of increasing 
the produced BAU by about 20%. 

Following Ref. [2H], in our numerical computations we only include processes in the 
thermal bath mediated by neutrino and top Yukawa couplings. This means that we 
neglect AL = 1 scatterings involving gauge bosons [301 120] , thermal corrections [20J and 
possible effects from spectator processes [3T] . but that we do take into account corrections 
from renormalisation group (RG) running between the electroweak scale and [El E2] 
(for which we use the software package REAP [33]). Regarding the pole mass of the top 
quark, we take the value mf ole = 174.2 GeV [HI]. We stress that the uncertainties in this 
value can have a strong influence on the RG evolution of the relevant parameters (namely 
the neutrino masses). Thus, the latter can provide a significant source of uncertainties 
in the BAU estimates. The renormalised value of the top Yukawa coupling (at energies 
J7ijvi), also directly affects the strength of the AL = 1 scatterings. Furthermore, let us 
notice that we also neglect AL = 2 scatterings, which is a good approximation as long 
as \A aa \m 1)a /m* > 10 x (m^/10 14 GeV) [32]. 

Finally, our estimates of the produced BAU are based on a set of coupled Boltzmann 
equations, which only provide a classical approximation to the Kadanoff-Baym equations. 
Quantum effects for thermalization have been ignored in our analysis |35j. Other approx- 
imations have led to the present simplified form of the Boltzmann equations: for instance, 
it is usually assumed that elastic scattering rates are fast and that the phase space den- 
sities for both fermions and scalars can be approximated as f(Ei,T) = {ni/n^)e~ E ^ T ', 
where = ^TmfK2(mi/T), with being the number of degrees of freedom. Accord- 
ingly, we also use the equilibrium number densities in this so-called Boltzmann approxi- 
mation. 

All the above mentioned approximations (as well as others we have not explicitly 
discussed) will ultimately translate in potentially relevant theoretical uncertainties when 
estimating the value of the BAU. Thus, it is important to bare in mind that one may be 
either over- or under-estimating ris/n^, so that caution should be exerted when deciding 
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upon the BAU viability of a given SUSY seesaw scenario. 

3.2 Constraints on the reheat temperature and the gravitino 
problem 

The predictions for the baryon asymmetry arising from a given seesaw scenario can 
be severely compromised due to constraints on Trh- One important example of such 
constraints emerges in locally supersymmetric theories, and stems from the decays of 
thermally produced gravitinos [3B1 E7J ■ In this class of SUSY models, and assuming a 
low-energy MSSM with R-parity conservation, either the gravitino is the lightest super- 
symmetric particle (LSP) and is thus stable, or else it will ultimately decay into the LSP. 
Two generic problems arise from these decays, and are the following. 

Firstly, gravitinos can decay late, after the Big Bang nucleosynthesis (BBN) epoch, 
and potentially spoil the success of BBN [3H1 EE]- This leads to upper bounds on the 
reheat temperature which depend on the specific supersymmetric model as well as on the 
mass of the gravitino [37] . For gravitino masses in the TeV region, the gravitino BBN 
problem practically precludes the possibility of thermal leptogenesis. However, with a 
heavy gravitino (roughly above 100 TeV), the BBN problems can be nearly avoided. If 
one considers the gravitino mass as a free parameter, one can safely avoid the latter 
constraints for any given reheat temperature. 

Secondly, the decay of a gravitino produces one LSP, which has an impact on the 
relic density of the latter. The number of thermally produced gravitinos increases with 
the reheat temperature, and we can estimate the contribution to the dark matter (DM) 
relic density arising from non-thermally produced LSPs via gravitino decay (for heavy 
gravitinos) as [38j E7] 

which depends on the LSP mass, mLSP, as well as on Trh- Taking the relic density 
bound f2LSp~ th ^ 2 < ^dm^ 2 < 0.13 from the Wilkinson Microwave Anisotropy Probe 
(WMAP) P, we are led to an upper bound on the reheat temperature of 

T RH < ,4xlO"Gev(^). (19) 

V. ™LSP / 

For masses of the LSP (assuming this to be the lightest neutralino) in the range 100 
GeV — 150 GeV we obtain an estimated upper bound on the reheat temperature of ap- 
proximately Trh < 2 x 10 10 GeV. Naturally, heavier LSP masses lead to more severe 
bounds on Trh- 
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There are other frameworks where, although still viable, thermal leptogenesis is signif- 
icantly constrained by bounds on T RH - This can occur for scenarios with stable gravitinos, 
i.e. where gravitinos are the LSP. In many cases the bounds on the reheat temperature 
strongly depend on the model under consideration, for instance on the properties of the 
next-to-lightest supersymmetric particle (NLSP). For example, recent studies of models 
with gravitino LSP have found the following bounds for the Constrained Minimal 
Supersymmetric Standard Model (CMSSM), 



while a scenario with gaugino-mediated supersymmetry breaking and sneutrino NLSP 
(stau NLSP) implies [ffi] 



In the subsequent numerical analysis, and as examples, we will take into account the 
following bounds, Trh < 2 x 10 10 GeV and Trh < 5 x 10 9 GeV, respectively represen- 
tative of a mildly and a strongly constrained case for thermal leptogenesis. 

4 Constraints on the seesaw parameters 

After having gone through the most relevant aspects of the computation of the BAU, let 
us now proceed to discuss the constraints on the several seesaw parameters which arise 
from assuming that the baryon asymmetry is generated by flavour-dependent thermal 
leptogenesis. 

The parameters of the Umns matrix, as well as the two mass squared differences, 
are presently constrained by neutrino oscillation experiments. From them we know that 
s 2 3 ~ C23 Rj a/1/2, su ~ a/1/3 (and cu ~ a/2/3), while #i 3 is only bounded from above, 
#i3 < 11.5° (at 3cr confidence level) [H]. Regarding the Dirac and Majorana phases (<5, 
if i and (£2), n o experimental data is yet available. For the case of hierarchical neutrinos, 
we have that m U2 « a/ Am^ « 0.01 eV and m V3 ~ a/ Am^ tm ~ 0.05 eV. On the other 
hand, parameters like the heavy neutrino masses m^, and the mixing angles involving 
the heavy Majorana neutrinos (i.e. the complex .R-matrix angles 9{), are experimentally 
unreachable. 

The main focus of this work is to address the constraints on mixing angles and CP 
violating phases of both the light and heavy neutrino sectors. However, and especially 
when reheating effects are taken into account, interesting bounds for the heavy and light 
neutrino masses can also be derived. We begin our discussion by briefly re-analysing the 
latter constraints for the case of the MSSM, considering flavour- dependent effects. 



7r H <4 x 10 9 GeV, 



(20) 



T RH < 7 x 10 9 GeV (T RH < 3 x 10 9 GeV) . 



(21) 
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4.1 Heavy and light neutrino masses 



Let us start with general considerations regarding bounds on the light and heavy neu- 
trino masses from thermal leptogenesis in the MSSM, when flavour effects are taken into 
account. From Eqs. (TI4"1 ITTI) . it is clear that within our framework (and approximations), 
the masses mjv 2 and mjv 3 are not constrained by thermal leptogenesis. 

In flavour- dependent leptogenesis, the decay asymmetries are constrained bj@ p2] 



£li „ < £2* m „ ( , ( ™*) . (22) 
87rfo \ mi J V m i / 



In the above equation, and for hierarchical light neutrino masses, e™ ax = 3 m U3 /(8iT i>|) 
is the maximally possible value, both in the flavour-independent approximation and in 
the flavour- dependent treatment. Regarding washout, in the type-I seesaw, the flavour- 
independent washout parameter satisfies 



mi > m Ul , (23) 

whereas the flavour-dependent washout parameters are generically not constrained. 
In the flavour-independent approximation, Eq. (I2"3l leads to a dramatically more restric- 
tive bound on E\ = ^2 a £i, a [21] . and finally even to a bound on the neutrino mass 
scale [32]. We also note that for quasi-degenerate light neutrino masses, an optimal 
washout parameteiEl rh\ >a ~ m* is possible, but it however implies that the decay asym- 
metries are suppressed by at least a factor yjm* /m Ul when compared to the optimal value 
£ max ^ c £ (J22]))J§ Concerning the decay asymmetries, we will see that the suppression 
factor yjrhi ta /rhi has further interesting implications also for the case of hierarchical 
neutrino masses. 

Using the upper bound on the decay asymmetry of Eq. ([221 for the case of hierarchical 
light neutrinos, and assuming an optimal efficiency, it is possible to estimate the baryon 
asymmetry. Even without reheating, the comparison of the estimated value with the 
observed BAU by WMAP [T], 

n B /n^ w (6.10 ± 0.21) x 10" 10 , or Y B w (0.87 ± 0.03) x 10" 10 , (24) 



4 For simplicity, we present the discussion for normal mass ordering. For inverse ordering, m Vz is 
replaced by m„ 2 . 

5 From here on, and regarding analytical discussions, we will assume |v4 Qa |mi !Q « fni, a , thus neglecting 
the 0(1) quantity A aa . The latter is included in the numerical computations. 

6 For instance, in the type-II seesaw, where an additional direct mass term for the light neutrinos from 
SU(2)L-triplets is present, this suppression can be avoided and the maximal decay asymmetry 3 "j Vl m U3 
(for normal mass ordering) can be realised for quasi-degenerate neutrinos [13] . It is easy to see that this 
holds true also in the flavour-dependent case. For instance, in "type-II upgraded" seesaw models (see 
e.g. [B]), this bound can be nearly saturated easily. 
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Figure 2: Isosurfaces corresponding to the maximally possible BAU as function of mjvi and Trh- From 
small to large rnjvi — Trh regions, the associated BAU ranges are: jib/w 7 € [10~ 10 , 3 x 10 -10 ], ns/n 7 S 
[3xl(T 10 , 5.9xlCT 10 ], n B /n 7 e [5.9xlCT 10 , 6.3xlCT 10 ], n B /n 7 e [6.3xlCT 10 , 1CT 9 ] andn B /n 7 > 1(T 9 . 
The corresponding colour code is grey, green, dark blue (WMAP), light blue and red, respectively. 

allows to obtain a lower bound on the mass of the lightest right-handed neutrino [21J. 

Clearly, the maximal BAU that can be generated depends on both and Trh- The 
combined constraints on these quantities are shown in Fig. [21 Leading to the latter, we 
have considered a normal hierarchical spectrum of light neutrinos. We have also assumed 
a maximal decay asymmetry as in Eq. (122 1) and an optimal efficiency for a given tunx/T^r. 
From Fig. [2l let us finally point out that in order to obtain BAU compatible with the 
WMAP range (represented in dark blue), the minimal values for the reheat temperature 
and for m Nl are T^ n « 1.9 x 10 9 GeV and m™ in « 1.5 x 10 9 GeV. Moreover, in the 
presence of an upper bound on the reheat temperature, there is also an upper bound on 
mjvi, stemming from the dramatic loss of efficiency occurring when ^> Trh, as it 
was shown in Fig. [TJ For instance, T RH ^2x 10 10 GeV imposes m™^ x 1.4 x 10 11 GeV, 
while T RH < 5 x 10 9 GeV yields m™^ x rs 1.9 x 10 10 GeV, leading to viability windows for 
the mass of the lightest right-handed neutrino. 

In the following analysis, whenever we present numerical examples regarding con- 
straints on the seesaw parameters, we will vary m^, and select the value for which the 
produced baryon asymmetry is maximal. Typically, this corresponds to around (or 
slightly above) Trh, as illustrated in Fig. [31 In order to reduce the BAU associated with a 
given choice of parameters, one can simply vary m^ l . Lowering such that <C m* 
leads to a regime where n^/n-y decreases with decreasing m^. On the other hand, it is 
also possible to increase m Nl , taking values m Nl ^> T RH , since then the strong washout 
leads to an important reduction in the BAU. 

Finally, and before concluding this discussion, let us comment on the bounds for the 
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Figure 3: Value of the ratio m^/Tjin which yields the maximally possible baryon asymmetry as a 
function of |v4 QQ |mi, Q . Solid, dotted, and dashed lines correspond to fhi/(\A aa \ fhi >a ) = 2,5 and 100, 
respectively. 



masses of the light neutrinos. At present, the absolute neutrino mass scale m Vl (for normal 
mass ordering) is only experimentally constrained by Tritium /3-decay, Gv/3/3-decay and 
cosmology, to be roughly below 0.5 eV pES] . In general, flavour-dependent leptogenesis will 
not provide any additional constraints on m Ul . As mentioned before, in the present study 
we focus on hierarchical light neutrinos. In this case, it is nevertheless interesting to point 
out that with a bound Trh < 5 x 10 9 GeV or Trh < 2 x 10 10 GeV, we have numerically 
verified that increasing the neutrino mass scale towards a quasi-degenerate light neutrino 
mass spectrum leads to a reduction of the BAU-allowed regions of the parameter space. 
This is essentially due to two reasons. Firstly, although there is no bound on m l Q! , 
its typical values are of the order of m ul ,m U2 and m U3 . Therefore, only for strongly 
hierarchical light neutrino masses can m ui set the right scale for an optimal washout 
parameter. Secondly, slightly increasing m Ul towards a quasi-degenerate spectrum does 
not significantly enhance e™ 8 *, but does reduce each of the decay asymmetries due to a 
factor of ^fhi a /fhi m \J m* jrh\ < \Jm* /m Ul , enforcing optimal washout m liQ « m* 
(using Eq. f[2"3]) ). This is illustrated in Fig. HI where we display the bound on the decay 
asymmetry ei jCt , normalised to the maximal decay asymmetry e max '° = \J Am^ tm 

for a (normal) hierarchical mass spectrum of light neutrinos. The washout parameter 
m ljCr has been fixed to m* = sin 2 /? x 1.58 x 10~ 3 (close to its optimal value). 



4.2 Mixing angles and CP phases 

Here we will discuss how the constraints arising from the reheat temperature can affect 
the washout and efficiency factors, and in turn favour/disfavour choices for the mixing 
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Figure 4: Bound on the decay asymmetry e\^ a as a function olm Ul . The washout parameter |-A QQ ,|mi iQ , 
is fixed to m* = sin 2 /3 x 1.58 x 1CP 3 (close to optimal), and the asymmetry is normalised to its maximal 
value e m "'°, obtained for a (normal) hierarchical mass spectrum of light neutrinos with m Ul = 0. The 
bound of the decay asymmetry in the type-I seesaw (solid line) decreases with increasing neutrino mass 
scale and remains below the bound for hierarchical masses for the depicted values of m Ul < 0.4 eV. For 
comparison, the general bound, which can be saturated in the type-II seesaw scenario, is also displayed 
(dotted line). The latter increases with increasing neutrino mass scale. We have used tan/3 = 30 and 
m Nl = 10 10 GeV. 



angles and CP violating phases. We also analyse the impact of the latter constraints 
regarding the flavour- dependent CP asymmetries, and investigate some illustrative limits 
regarding the i?-matrix angles, 

4.2.1 Favoured parameter regions with optimal washout 

With strong constraints on the reheat temperature like Trh < 5 x 10 9 GeV or Trh < 
2 x 10 10 GeV (motivated by Eqs. (jl9H21|) ). and a lower bound on m Nl of about 10 9 GeV, 
only a rather small window of values remains allowed (c.f. Fig. [2D- This means 
that in order to produce enough baryon asymmetry, at least one of the efficiency factors 
r] a = T)(A aa mi ia , mi, TRH/rajVi) should be close to optimal (i.e., mi )Q must not differ much 
from m*) and, in addition, the corresponding decay asymmetries should approach e™ ax . 

An optimal efficiency r\ can be achieved for rhi t0l ~ m* = sin 2 j3 x 1.58 x 10 -3 eV. Note 
that m* is much smaller than m U3 and, even though m U2 < m U3 , m* is still significantly 
smaller than m U2 . For hierarchical light neutrino masses, apart from having m vi m U2 < 
m U3 , m Vx is unconstrained and can vary in a large range. For instance, it is possible to 

7 We would like to note that our model- independent results, presented in terms of the i?-matrix 
parameterisation, can be understood as constraints on whole equivalence classes of neutrino mass models 
(see e.g. US]). 
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choose m vi close to m*, such that m l Q m m Vl would ensure an optimal value for the 
washout parameter. In the following numerical analysis, we will set m ui = 10~ 3 eV ~ m* 
as an illustrative example. 

To begin the analysis, let us start by writing the washout parameters explicitly in 
terms of the seesaw parameters, using the i?-matrix parameterisation. One has 

mi, e = |e2 vi c 2 c 3 c 12 c 13 y/m^ + s 3 (e lS s 1 s 13 y/m^- e^ 2 c x s i2 c i3 

'X i I 2 

+s 2 c 3 (-e l ci S13 - e2^ 2 si S12 ci 3 v /m^") | , (25) 

fhl,» = | ~ e^ 1 C2C 3 S 1 2C2 3 ^/^ + S 3 (siS 23 y/^ - e^ 2 C 1 C 1 2C2 3 ^/^) 

+s 2 c 3 (-c 1 s 23 y/^- e^ 2 s 1 c 12 c 23 y/m^) + 0(8 13 )\ , (26) 
rhi,r = \e^ x c 2 c 3 s 12 s 23 y/m^ + s 3 (s 1 c 23 y/m^ + e^ 2 a c 12 s 23 y/m^) 

i 1 2 

S2C 3 (-C 1 C 23 y/ ; ^+ 6^ 2 S 1 C 12 S2 3 y/m^) + 0{8 V3 )\ . (27) 

For fhi /j,, fhi jT , only the zeroth order terms in an expansion in 8\ 3 are shown. We note 
that, when compared to 82 and 8 3 , 61 only plays a minor role for leptogenesis. This can 
be clarified by considering the second and third terms on the right side of the equations, 
which contain the two potentially large contributions associated with y/m U2 and y/m U3 . As 
can be seen, a real 8\ mainly rotates terms proportional to y/vn V2 into terms proportional 
to y/m u:j , and thus it is not introducing any new features. For instance, with 81 = n/2, 
the roles of 8 2 and #3 are simply interchanged: a 8 2 rotation now generates contributions 
to m lQ , of 0{m V2 ) and a 8 3 rotation induces contributions to rn,i )Q of 0(m U3 ). Large 
imaginary parts of 8\ typically lead to fh\ >a > m Ul « m*, and thus are not useful in 
achieving optimal efficiencies. Although it is straightforward to generalise the discussion 
to include arbitrary 8\, we will simplify the analysis by setting 8\ = in what follows. In 
the limit of 8\ = 0, we find 

mi, e ~ |e5 Vl C2C3Ci3Ci2 v /^-e^ 2 S3Si2Ci3 1 /^-e j5 S2C 3 Si3 1 /r^| 2 , (28) 
1 * * 1 2 

Wi i(U w I - e^ 1 C2C3S12C23 y/rri^- e^ 2 s 3 c 12 c 23 y/m^ - s 2 c 3 s 23 y/m^\ , (29) 

I i i I 2 

mi, T ~ I e 5 Vl c 2 c 3 sia s 23 v/^v + e^ 2 s 3 c 12 s 23 y/m^ - s 2 c 3 c 23 y/m^\ , (30) 

where, as in Eqs. ( 12"B"1 127]) . only the zeroth order terms in an expansion in # 13 are shown 
for mi i(1 ,m^ T . In general, for an optimal efficiency, it is crucial to avoid contributions to 
fhi )Ct of order 0(m U2 ) or 0(m U3 ). In the flavour-independent approximation, achieving 
an optimal washout (given by the sum rh\ = ^ a wzi, a ) necessarily required that both 
\s2\ and I S3 1 were small |19j . As can be seen from Eqs. (I28H30I) . when flavour effects are 
included, washout considerations still favour parameter regions with small | S2 1 and | S3 1 . 
Nevertheless, in the flavour- dependent treatment, the individual washout parameters rh\ ta 
are typically smaller than their sum, mi, and can significantly differ from each other. It 
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is therefore pertinent to re-investigate whether other regions of parameter space may also 
allow for optimal washout. 

Let us first consider the contribution to m 1)T (or similarly to rhi ^) proportional to m U2 , 
due to non-zero values of 63. From Eq. (l30l) we find that this is given by c\ 2 s\ m V2 ~ 
I s§ m V2 . In addition, the quantity which enters the efficiency factor is the product 
A TT mi jT , with A TT « —0.6. Combining these two effects reduces the ^3-induced washout 
by a factor of about 1/5, and even for s 3 ~ 0{1) optimal washout |A aQ |m liCf « m* 
is still possible to obtain. On the other hand, the contribution to m l r (or similarly 
to fh 1)fJ ) proportional to m U3 due to non-zero 9 2 is given by c^ z s\c\m V3 ~ \s\c^m vz . 
Again, the quantity which enters the efficiency factor is the product |74 rr |mi jT . Although 
combining these two effects reduces the #2-induced washout by a factor of about 0.3, 
optimal washout cannot be achieved for large | ^ 2 U if ^3 i s small. The only exception 
occurs for small |c 3 |, which can suppress the large contribution of 0(m U3 ) to m l r ,m 1)M 
and still allow for large | s 2 1 - However, as we will see in the next subsection, the decay 
asymmetries Ei j(X are somewhat suppressed in this case. 

Another difference between the flavour-independent approximation and the correct 
flavour- dependent treatment becomes apparent when we consider fhx >e . In contrast to 
mi )T and m^, large values of S2 only induce a washout parameter mi fi of sf 3 
In other words, we can be in the optimal washout regime m 1)e « m* also for larger values 
of s 2 , depending on the size of 9i 3 . Nevertheless, in this case the decay asymmetry e l e is 
also suppressed when compared to the optimal regime. 

Let us point out that another possible way to obtain optimal rhi^ and m l r of order 
m* would be to align s 2 and S3 such that their contributions to one of the washout 
parameters (respectively proportional to m V3 and m U2 ) would nearly cancel each other. 
In this case one could obtain, for example, rh l fl w m*. However, one can verify that 
the third washout parameter m l T (and thus rhi) would still be 0(m u:j ) or 0(m U2 ), which 
would imply a suppression of ei tft by a factor m*/fhi, in comparison to e™ ax (c.f. Eq. ([22]) ). 

4.2.2 Flavour-Dependent Decay Asymmetries 

From the analysis of the flavour-dependent washout parameters m 1]tt , it has become 
apparent that there are several regions of the seesaw parameter space with appealing 
prospects for leptogenesis. They can be summarised as follows. Generically favoured is 
the region of small \s2\ and \ss\, where the washout parameters rhi ja are 0(m Ul ), and 
receive only small contributions proportional to m V2 and m v , A . In the flavour-independent 
approximation, this was the only region favoured by washout [19J. In contrast, in flavour- 
dependent leptogenesis, washout (in all flavours) remains optimal also for larger values 
of I S3 1 < 1. With small |c 3 | (implying large \s^\ which is no longer disfavoured), large | s 2 1 
becomes compatible with optimal washout, and a whole new region, washout-favoured, 
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has emerged. An additional interesting effect is that washout in the e-ffavour, governed 
by fhi >e , can also be optimal for large \s 2 \, due to the smallness of #13 . 

In addition to an optimal and given the tight constraints on the reheat temper- 
ature, it is also desirable to have decay asymmetries close to the optimal value ef ax . Let 
us now address the latter issue more thoroughly. First, we note that for R — 1, and more 
generally, for Qi = 0, ~ (mod ir), the decay asymmetries exactly vanish. A deviation from 
the latter values of the form 6 2 = 3 = but 6i ^ 0, also leads to zero values. To clarify 
the analysis, let us explicitly write the decay asymmetries in terms of the seesaw parame- 
ters. For simplicity, and as illustrative examples, we first consider the dependence on 8 2 , 
with 9 3 = Q\ = 0, and then study the ^-dependence, setting 9 2 = 9\ = 0. To discuss the 
parameter regions with optimal washout where 9 2 and #3 are both large, we then turn 
to the decay asymmetries with nonzero 9 2 and #3. For simplicity, in the latter case we 
will again present the formulae with 6\ = 0. The discussion can be easily generalised to 
arbitrary 8i values, using the general expressions for the decay asymmetries as given in 
Eqs. (fT2l [I~4"l) . We will also illustrate the effects of nonzero 9\ via numerical examples in 
Section 14.31 

For #3 = 61 = 0, one obtains 
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(31) 



(32) 



(33) 



■ 

max 



where the dots indicate terms 0{mt x y/m V3 ) and ©(m 2 ^), which are sub-leading for hier- 
archical light neutrinos. Eqs. (152] 1331) show that, provided that 6 2 has an imaginary part 
the ^2-induced decay asymmetries and e ljT can be close to the optimal value e" 
defined in Eq. (|2*2*|) . Regarding the sub- leading terms in Eqs. (I3"2"| l3"3l . decay asymmetries 
can also emerge if the Majorana CP phase tpi is nonzero, albeit suppressed by a factor 
0(y/m Ul /m V3 ). On the other hand, for hierarchical light neutrinos, the leading contri- 
butions to £i )f . are suppressed by s 2 3 and S13 when compared to and e 1)T . This is in 
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agreement with Eq. f[2"2"]) . which states that for large |s2| (and thus rh\ ~ 0(m U3 )) the 
decay asymmetry is suppressed by a factor \fm^Jm{ s r3 , compared to the maximal 
value, e™ ax . 

Let us now consider the case where 62 = Q\ = 0. The flavour- dependent decay 
asymmetries are as follows: 
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(36) 



In the above, the dots denote sub-leading terms 0{ml x y/m Ua ) , 0{m^^). Compared to 
ef ax , the decay asymmetries induced by 6 3 are suppressed by a factor of 0(m U2 /m U3 ) « 
1/5 with respect to those induced by 6* 2 . As can be seen from Eqs. ( I3"4"ll3"6"l) . either a non- 
zero imaginary part of 6 3 or non-zero phase 5 or the difference ifi — <f2, can in principle 
provide the required CP violation for leptogenesis. 

In our discussion of the seesaw parameters with optimal washout, we have encountered 
a new region where |s 3 | and |s2| were large, but with small |c 3 |. In the limit of vanishing 
6\, and keeping only the leading terms for simplicity, we obtain 
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In the above, the dots denote sub-leading terms 0{mi 1 ^m U3 ), 0(mj} 2 ^m U3 ), C^m^J and 
0{m 2 U2 ). Eqs. (I37H39P show that for large |s 3 |, the asymmetries are suppressed by |c 3 |/|s 3 | 
when compared to the optimal value e™ ax defined in Eq. ( !22l) . Although it is possible to 
achieve optimal washout in the region with large |s 3 | and |s2| if |c 3 | is small, the smallness 
of |c 3 | in turn leads to somewhat suppressed decay asymmetries in all flavours. 

To conclude this discussion, let us stress that as seen in the above considered cases (c.f. 
Eqs. ( 13TH39]) ) . and contrary to what occurred in the flavour- independent approximation, 
even for a real i?-matrix, one can indeed obtain non-vanishing values for the baryon 
asymmetry. Notice however that the contributions of the £7mns phases appear suppressed 
by ratios of the light neutrino masses (^m Ui /m Uj , i < j) and/or by #i 3 . Also notice that 
other new regions of the parameter space where \s2 \ and |s 3 | can be large (with complex R) 
have good prospects for leptogenesis, but lead to decay asymmetries which are suppressed 
when compared to the favoured region of small | ^ 2 1 and |s 3 |. 

In the following subsection, we will discuss whether or not the promising regions here 
identified still offer viable BAU scenarios when bounds on Trh (namely Trh < 2 x 10 10 
GeV and Trh < 5 x 10 9 GeV) are taken into account. 

4.3 Numerical examples 

The analysis of the previous subsections, based on the analytical expressions for the 
flavour dependent washout parameters m l cr and decay asymmetries £i )Q , has revealed in- 
teresting differences between the constraints on the seesaw parameter space arising from 
leptogenesis in the flavour-independent approximation and the correct flavour- dependent 
treatment. In what follows, we will present some numerical examples which illustrate 
these new constraints, taking into account bounds from Trh- In all the following exam- 
ples, we will always consider m Vl = 10~ 3 eV. 

Before beginning, let us notice that the figures here displayed show ranges of the 
maximal attainable BAU. As done before, a scan is performed over m^ 11 and its value is 
determined as to obtain maximal n^jn^. Additionally, it is important to stress that, in 
agreement with the discussion of Section l3.1.3l there may still exist significant theoretical 
uncertainties in the estimates of the produced baryon asymmetry. As previously men- 
tioned, the effect of these uncertainties is hard to quantify, and can lead to both over- 
and under-estimations of the BAU. An educated guess of these theoretical uncertainties 
would suggest that one should allow for as much as a factor 2 (or even 5) between the 
real and the estimated values. Thus, when evaluating the BAU viability of the seesaw 
parameter space, we will also be showing regions where the produced baryon-to-photon 
ratio lies outside the WMAP observed range, but is still larger than 10~ 10 . In particular, 
we allow for a factor 2 (5) arising from theoretical uncertainties, and display the corre- 
sponding regions n^jn^ G [3 x 10~ 10 , 5.9 x 10~ 10 ] (n B /n^ G [10~ 10 , 3 x 10~ 10 ]) in green 



22 



(grey). If our computations were exact, the shown region with n^/n^ > 5.9 x 10~ 10 could 
be compatible with the observed baryon asymmetry (notice that values larger than the 
WMAP range can be easily accommodated by varying m^). 

Even in the absence of considering the new CP violation sources arising from the C/mns 
matrix, the flavour- dependent computation gives rise to interesting new constraints on 
the seesaw parameter space. Thus, we first examine a conservative scenario with CP 
violation exclusively stemming from complex i?-matrix angles, taking into account the 
bounds for the reheat temperature. 

Figure E] illustrates the Ke(9 2 ) — lm(9 2 ) and Re(6 , 3 ) — Im(# 3 ) regions compatible with 
successful thermal leptogenesis in the presence of bounds T RH < 2 x 10 10 GeV and 
2~rh < 5 x 10 9 GeV. In this case 9i3,5,(pi,(p 2 and 9\ have been set to zero. On the left 
(right) panels, 9 3 = (0 2 = 0). The examples with T RH < 2 x 10 10 GeV (i.e. Fig. E^a) 
and Fig. [3(b)) update the analysis of Ref. [19], which had been performed in the flavour- 
independent approximation. In the present flavour-dependent computation, we find that 
the BAU arising from complex # 3 is somewhat larger and new regions, which are indeed 
compatible with the WMAP range, have now emerged. It is worth stressing that in this 
case of complex R-mathx angles, the favoured regions still correspond to small values 
of 9 2 and #3. Considering a stronger bound on Trh, namely Trh < 5 x 10 9 GeV, we 
notice that there are still regions in the ~Re{9 2 ) — \m.{9 2 ) plane compatible with WMAP 
observations. When the latter bound on T RH is applied, we verify that for complex 9 3 
values it is no longer possible to saturate the WMAP preferred range. Nevertheless, 
regions where ub/u^ G [10~ 10 , 3 x 10~ 10 ] can still be found (viable if one allowed for a 
factor 5 uncertainty in the computation). In any case, it is manifest that for this stricter 
Trh bound, the preferred source of CP violation for leptogenesis is 9 2 . In both cases, 
the observed differences between the present and the previous analyses (|19j) originate 
from taking into account flavour effects in the Boltzmann equations. In the 9 2 case, the 
differences are less apparent, essentially due to the fact that fhi ie = and fhx ttM ~ fhi jT . 
However, important effects can be observed for the Re(# 3 )— Im(#3) plane, since in this case 
both the decay asymmetries and washout parameters differ for each individual flavour. 
In particular, this leads to deformations of the allowed regions when compared to those 
presented in Ref. [T9"] . 

In Fig. [U] we illustrate the effects of non-zero 9%. Taking arg(6* 2 ) = arg(6* 3 ) = 7r/4, and 
again 9±s, 8, ipi, if 2 = 0°, we now display the regions of the Re(6*2) — Re(6' 3 ) parameter 
space compatible with thermal leptogenesis in the presence of a bound T R h < 5 x 10 9 
GeV. First, let us point out that Fig. [6]^a) corresponds to a variation of Figs. [5](c) and 
(d), but for fixed values of arg(6* 2 ) and arg(0 3 ). When compared to Fig. 0(a), Fig. Wh) 
shows the effect of 9\ = 7r/4, which is mainly a rotation (and a slight deformation) of 
the allowed region. Figure EJ^c) and Fig. [6]^d) illustrate that an imaginary part of #1, in 
addition to introducing an additional source of CP violation, leads to a reduction of the 
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Figure 5: Regions of Re(6>2) — Im(&a) and Re(#3) — Im(03) parameter spaces compatible with successful 
thermal leptogenesis in the presence of bounds on the reheat temperature: Trh < 2 x 10 10 GeV (upper 
panels) and Trh < 5 x 10 9 GeV (lower panels). From out- to inner- most regions, the associated ranges of 
maximally possible baryon asymmetry are: ne/n 7 G [10~ 10 , 3 x 10 -10 ], ne/n 7 G [3 x 10~ 10 , 5.9 x 10~ 10 ], 
n B /n 7 G [5.9 x 1CT 10 , 6.3 x 1CT 10 ], n B /n 7 G [6.3 x 10" 10 , 10" 9 ] and n B /n 7 > 1(T 9 . The corresponding 
colour code is grey, green, dark blue (WMAP), light blue and red. 
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compatible parameter space. As argued in Section 14.2.14 this effect can be explained by 
a stronger washout due to an enhancement of the parameters rhi )0l . We emphasise that 
for complex Q{ (but vanishing low-energy CP phases), even in the presence of non-zero 
values of 6%, small values of |s2| and \s%\ are still favoured by thermal leptogenesis. 

In Fig. [5] we have separately considered the effects of each of the i?-matrix angles 
#2,3, while in Fig. [6] we analysed the impact of non-vanish 8\ upon the Re(#2) — Re(#3) 
parameter space, assuming sizable arguments for both 8 2 and #3. However, for quite small 
values of the arguments, and when flavour effects are taken into account, new interesting 
regions of the ~Re(8 2 ) — Re(6 , 3 ) parameter space can also arise. This is shown in Fig. [7J, 
where we now display the regions of the Re(#2) — Re(^) parameter space compatible with 
successful thermal leptogenesis, assuming a small value for both arguments, namely it/ 16. 
As can be seen, not only do we encounter large values of BAU associated with small values 
of |s 2 | and |s 3 |, but new extensive regions, with larger values of 8 2 and 9 3 , are now present. 
The origin of these new regions exhibiting a sizable BAU can be easily understood (in 
the limit of vanishing 81) from the analytical considerations of Sections 14.2.11 and 14.2.21 
On the one hand, Eqs. (12811301) show that optimal washout is possible in two cases: for 
large S3, or then for large S3 and s%, provided that C3 is small. On the other hand, from 
Eqs. (I37H39P we have seen that optimal decay asymmetries require contributions from 
non-zero 82, suppressed if C3 is small. The shape of the extensive regions in Fig. [7Ja) 
with sizable BAU reflects the balance between having a sufficiently small washout, while 
at the same time succeeding in obtaining an important decay asymmetry. As expected, 
taking into account stricter bounds on Trh leads to the disappearance of the WMAP 
compatible regions (c.f. Fig. [7J(b)). Nevertheless, regions where the BAU is close to the 
observed range still survive, both for small and large 82 and 83. For non- vanishing values 
of 6\, just like discussed regarding Fig. El one would observe a deformation of the regions 
displayed in Fig. [7J The analytical interpretation can be now obtained from Eqs. ( [2511271) . 
albeit in a less straightforward way. 

After having revisited leptogenesis scenarios where CP violation originated solely from 
the complex i?-matrix angles, let us now consider the effects of having CP violation 
arising from the £/mns phases. This is especially appealing given that, and contrary to 
the i?-matrix angles, parameters like S and #13 are likely to be observable in neutrino 
oscillation experiments. Additionally, and as pointed out in Refs. [T5l IT6] (although 
in the context of the SM), these are examples of scenarios where there is a maximal 
connection between leptogenesis and low-energy CP violation. We begin by addressing 
a scenario where the i?-matrix angles are real (but non-zero) and 5 is the only source of 
CP violation. For non-zero #13 and 5, Fig. [H] illustrates the emergence of new regions in 
the Re(# 2 ) — Re(^) parameter space, potentially compatible with thermal leptogenesis in 
the presence of a bound Trh < 2 x 10 10 GeV. In this example we have chosen 9i — tpx — 
if2 = 0°. Notice that in the flavour-independent approximation, leptogenesis would have 
been impossible for a real i?-matrix. From Fig. [HI we find that for 6* 13 = 11.5° (the largest 
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Figure 6: Regions of Re(# 2 ) — Re(#3) parameter space compatible with successful thermal leptogenesis 
in the presence of a bound Trh < 5 x 10 9 GeV. We display several values of complex 9\, choosing 
arg(# 2 ) = arg(^3) = 7r/4 and 0x3 — 8 — — — 0°- Colour code as in Fig. [5] 
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Figure 7: Regions of Re(# 2 ) — Re(# 3 ) compatible with successful thermal leptogenesis. We take B\ = 0, 
and consider arg(6> 2 ) = arg(6> 3 ) = tt/16. T rh < 2 x 10 10 GeV in (a), while T RH < 5 x 10 9 GeV in (b). In 
this example we have set #13 = S = <p\ = (f2 = 0°. Colour code as in Fig. [5] 



value experimentally allowed) and CP violating phase 5 close to 7r/2 (which maximises 
the decay asymmetry £i, e ), somewhat larger 82 values could now be marginally allowed. 
In any case, the largest values of the BAU are still associated with small 82 and 83. On the 
other hand, moving away from the present upper bound on # 13 , we find that the scenario 
is even more compromised. In fact, for 13 = 7.5°, only regions with BAU differing from 
the WMAP range by a factor 5 survive. For smaller #13 values (namely #13 < 5°), n^jn^ 
is already well below 1CT 10 . Likewise, considering stricter bounds on Trh would lead to 
the disappearance of all the shaded regions of Fig. 

In Fig. [9l we display the effect of non-zero 8\ on the allowed regions with non-zero #13 
and S. In the first example with 8\ = tt/4 and #13 = 11.5° (Fig. E^a)), we observe that in 
addition to a rotation of the allowed parameter space, the latter is somewhat enlarged. 

At this point, and from the examples so far considered, we are led to the conclusion 
that, unless # 13 is found to be close to its present upper bound, it is quite difficult to 
accommodate viable BAU scenarios relying on 5 as their only source of CP violation. For 
stricter bounds on the reheat temperature, the latter scenarios become increasingly more 
compromised. 

In addition to 5, there are other sources of CP violation arising from the [/mns matrix, 
namely the Majorana phases (pi, (p2- in principle, the latter could also provide the required 
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Figure 8: Regions of Re(6 l 2) — Re(6*3) parameter space compatible with successful thermal leptogenesis 
with a bound Trh < 2 x 10 10 GeV. In this example we have chosen real i?-matrix angles and 9\ = <pi = 
(fi2 = 0°. Colour code as in Fig. [5] 




Figure 9: Example displayed in Figs.[U[a) and (c), but with 6\ ^ 0°. Colour code as in Fig. [5] 
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Figure 10: Regions of Re(#2) ~ Re(#3) parameter space compatible with successful thermal leptogenesis 
with Trh < 2 x 10 10 GeV. In this example we have chosen real R, 0\ = and 6*13 = 5 = 0°. Colour code 
as in Fig. [5l 



CP violation for leptogenesis (see also Refs. [ISl HS1 US] ) • For completeness, in Fig. [TUlwe 
separately illustrate their role in generating a non-vanishing BAU. To do so, we assume 
a real i?-matrix, 5 = 0°, and take <pu2) = 0° on the left (right) panel. As seen from 
Fig. [10], when CP violation is exclusively arising from the Majorana phases it is indeed 
possible to obtain marginally compatible BAU values. Again, the most promising regions 
appear associated with small 9 2 and 6*3. We also observe that tpi and <p 2 lead to somewhat 
distinct regions on the Re(#2) — Re(#3) parameter space. In this example we have again 
considered a more relaxed bound for the reheat temperature, Trh < 2 x 10 10 GeV. 
As occurred for the cases investigated in Fig. El stronger bounds on Trh would imply 
that the generated BAU would also lie below 10~ 10 , so that the shaded regions of the 
Re(6 l 2 ) — Re(6 , 3 ) parameter space displayed in Fig. [TU] would disappear. 

Albeit it is pedagogical to consider the individual role of each Umns phase regarding 
BAU, in the most general case 5, ipx and ^2 can be simultaneously non-vanishing. In 
fact, and as shown in Fig. [TT], the Majorana phases can slightly improve the BAU allowed 
regions associated with 6 (and #13). We recall that, as discussed in relation with Fig. El 
values #13 < 5° failed to induce > 1CU 10 . Comparing Fig. CCD with Fig. Etc), we 

observe that for the choice of (pi 2 = 180° the regions where n^/n y 6 [3 x 10~ 10 , 5.9 x 
10~ 10 ] have become larger. Even though the WMAP range cannot be accounted for, it 
is nevertheless clear that the joint effect of the different C/mns phases translates in an 
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Figure 11: Same examples as in Fig.[]I)l but with 613 = 10° and S = 90°. 



improved scenario. 

Even though flavour-dependent thermal leptogenesis opens the possibility to generate 
the observed BAU exclusively from the £/mns CP violating phases, this may not be 
the most general nor the most successful scenario. The analysis of the present section 
lends a strong support to the latter statement. As we have found, it is quite difficult to 
encounter viable BAU scenarios associated with only low-energy CP violation. Moreover, 
if the given SUSY model implies a more stringent bound on the reheat temperature, BAU 
solely from C/mns phases becomes almost inviable. 

Recall that in the most general case, the R- matrix angles are also complex and that, as 
seen from Figs. 0, [U] and there are important regions in the 62-G3 parameter space where 
one can easily have compatibility with the WMAP range. In addition, i?-matrix phases 
allow for viable BAU even under a bound Trh < 5 x 10 9 GeV. Thus, it is important to 
investigate the simultaneous effect of all CP violating phases. In particular, one wonders 
to which extent the £7mns phases can affect the BAU predictions from the i?-matrix 
phases, and vice-versa. 

In Fig. [T2] we display the outcome of taking, in addition to the CP sources considered 
in Fig. [HJ non- vanishing values for the arguments of 62 and #3, namely 7r/16 (upper) and 
7r/4 (lower). The results are again depicted in the Re(6*2) — Re(6*3) plane, and we assume 
two bounds for the reheat temperature, Trh < 2 x 10 10 GeV (left) and Trh < 5 x 10 9 GeV 
(right). It is manifest from the comparison of Fig. [7] with Fig. [12] that the predictions 
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for the Re (#2) — Re (63) plane are hardly affected by considering non- vanishing values for 
the Umns phases. In the case of arg(6 l 2i 3) = 7r/16, and even for nearly maximal values of 
#13, one only observes a small distortion of the regions associated with small |s2| and | S3 1 , 
and a deformation of the regions associated with large | S2I and [S3 j (c.f. Fig. [T27b)). For 
larger arguments of 62 and 63 (7r/4), comparing Fig. ITWd) with Fig. EJ^a) implies that #13 
and S have had virtually no effect on the shape of the BAU compatible region, which is 
determined by the right-handed neutrino complex parameters 82 and #3. Finally, notice 
that for Trh < 5 x 10 9 GeV, it is possible to nearly reach the WMAP range for arg(#2,3) 
as small as 7r/16. Assuming larger values for the arguments of the i?-matrix angles allows 
to encounter wider regions where one still has compatibility with WMAP observations 
(panel (d)). However, this again favours the region of small [ S2 j and | s 3 1 . 

In a sense, Fig. [12] provides an illustrative summary of our analysis. Firstly, it con- 
firms the supposition that flavour effects in the Boltzmann equations indeed lead to the 
occurrence of new regions where BAU is viable (due to the already mentioned reduced 
washout). In addition, from its comparison to Fig. El it is also manifest that, in general, 
the leading role in BAU appears to be played by the i?-matrix complex angles, and not 
by the Umns phases (which in turn may compromise the possible bridges that could oth- 
erwise be established between low-energy CP violation and BAU). Moreover, Fig. [T2]is a 
clear example of the impact of the reheat temperature in severely constraining the SUSY 
seesaw parameters. The stricter the bounds on Trh, the more favoured are regions with 
small Re (6*2,3) and large arg(0 2 ,3)- 

5 Summary and conclusions 

In this study we have investigated the constraints on the SUSY seesaw parameter space 
arising from flavour-dependent thermal leptogenesis in the MSSM in the presence of 
upper bounds on the reheat temperature of the early Universe. In the temperature range 
here considered, both tau- and muon-flavours are in thermal equilibrium, so that the full 
flavour- dependence was taken into account. In order to calculate the efficiency factor for 
thermal leptogenesis, we have extended the flavour- dependent Boltzmann equations [TQl 
[TTl [T2"] . which were adapted to the MSSM case in [H], to include reheating (following 
Ref. [20J). Parameterising the solutions to the seesaw equation by means of a complex 
orthogonal matrix R [7], we have analysed which regions of the seesaw parameter space 
generically enable optimal efficiency and/or optimal decay asymmetries for leptogenesis. 

We have discussed several differences between the flavour-independent approximation 
and the correct flavour-dependent treatment of thermal leptogenesis. These are extensive, 
and, together with the bounds from the reheat temperature, lead to interesting new 
constraints on the SUSY seesaw parameter space. 
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Figure 12: Regions of Re((9 2 ) — Re^) parameter space compatible with successful thermal leptogenesis. 
We take 9 13 = 10°, S = 90°, and ip x = 180°. We consider arg(6> 2 ) = arg(6> 3 ) = tt/16 for panels (a) 
and (b), and arg(0 2 ) = arg(6> 3 ) = tt/4 for panels (c) and (d). T RH < 2 x 10 10 GeV ((a) and (c)) while 
Trh < 5 x 10 9 GeV ((b) and (d)). Colour code as in Fig. [SJ 
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Considerations on Trh give rise to the first constraints, in the sense that a dramatic 
drop in the efficiency takes place for ^> T RH (as much as three orders of magnitude 
for « 10 Trh). Since we have assumed that only the decays of the lightest right- 
handed neutrino were relevant for the lepton asymmetry, no bounds on the masses mjq 2 
and Tti were derived. Assuming an optimal regime for the decay asymmetry and for 
the efficiency, as well as that the inflaton only decays into MSSM particles (and not 
directly into right-handed (s)neutrinos), the requirement of a successful BAU leads to 
lower bounds on mjvi as well as on the reheat temperature. In particular, we have found 
m^ n k 1.5 x 10 9 GeV and Tjgf ~ 1.9 x 10 9 GeV, similar to the results obtained in the 
flavour- independent approximation [20]. On the other hand, in the presence of upper 
bounds on the reheat temperature (from dark matter relic abundance considerations), an 
upper bound on can also be inferred. In order to illustrate the impact of reheating, we 
have considered two examples for Trh, corresponding to mildly and strongly constrained 
scenarios: Trh ^2x 10 10 GeV, and Trh ^5x 10 9 GeV. Regarding the upper bound on 
m Nl} the latter bounds respectively yield m^ x « 1.4 x 10 11 GeV, and m™^ fa 1.9 x 10 10 
GeV. This leads to viability windows for the mass of the lightest right-handed neutrino. 

Regarding the light neutrino masses, namely m Ul , in general there is no upper bound 
from flavour- dependent thermal leptogenesis. Nevertheless, increasing m Ul within the 
present allowed experimental upper bounds (of roughly 0.5 eV), generically results in a 
reduced BAU. Furthermore, in the presence of strong Trh bounds, quasi-degenerate light 
neutrino masses (via the type-I seesaw mechanism) become disfavoured. 

As in the flavour-independent approximation, considerations on the washout param- 
eters generically favour the region of small | sin ^2 1 and | sin ^3 1 . However, and in clear 
contrast to the flavour-independent approximation, new regions with optimal flavour- 
dependent washout parameters have emerged, in association with large values of | sin #2 1 
and I sin #3 1 . In any case, in order to produce sufficient BAU in the presence of mild (or 
even strong) constraints on the reheat temperature, the decay asymmetries must also be 
close to their optimal values. Regarding the flavour- dependent decay asymmetries, in the 
general case of complex R matrix (but considering the limit of a real £7mns)> we gener- 
ically recovered the main results of the one-flavour computation (in the sense that the 
favoured regions still corresponded to small | sin# 2 | and | sin 6* 3 1 , albeit slightly enlarged). 
Another important result of our analysis concerns the effects of the reheat temperature, 
which were clearly manifest. In fact, taking stronger bounds on Trh leads to a significant 
reduction in the BAU allowed regions of the Re(6 l 2)-Im(6 l 2) and Re(6 l 3)-Im(6 l 3) parameter 
spaces (even to the disappearance of the WMAP compatible regions). In particular, for 
Trh < 5 x 10 9 GeV, we have seen that 9 3 cannot exclusively account for the observed 
WMAP results. 

In flavour-dependent leptogenesis, a potentially important role can also be played by 
the Uutsss phases. In principle, viable BAU scenarios could be obtained in the presence 
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of a CP-conserving .R-matrix, with the required amount of CP violation stemming either 
from the Dirac phase 5, or from the Majorana phases <pi and <p2- In the SM, this situation 
has been discussed in Refs. [T5| [T6] . However, the constraints on the seesaw parameters 
in the MSSM are expected to differ from the SM case, since for the temperatures (and 
values of tan f3) under consideration, all flavours are separately treated in the MSSM 
Boltzmann equations, whereas only the r-flavour is separately considered in the SM case. 

Exclusively relying on the phase 5 and on # 13 (under the standard parameterisation 
of the £/mns matrix) is a phenomenologically challenging choice, since these are the most 
likely (yet) unknown £/mns parameters to be experimentally measured. However, we 
have verified that even with 5 = tt/2, for values of #13 < 11.5° (the present experimental 
limit) the obtained values of the baryon asymmetry are only marginally compatible with 
observation (when large theoretical uncertainties are allowed for). By themselves, and 
even in the limit #13 = 0°, both Majorana phases, ipx and ip 2, could in principle account 
for BAU. However, and similar to 5, only marginal consistency with observations can 
be reached. In both cases the impact of the reheating temperature becomes manifest, 
since lower values of T RR can render these scenarios inviable. We also note that in the 
presence of small arg(#j), the BAU generated from CP violation in the right-handed 
sector dominates over the contributions from low-energy phases. Thus, the sensitivity to 
the £/mns CP violating phases is lost. In the limit of very strict bounds on the reheat 
temperature, one is thus compelled to take into account complex R as an additional 
source of CP violation. 

In summary, we have investigated which regions of the SUSY seesaw parameter space 
are favoured by flavour-dependent thermal leptogenesis, when bounds on the reheat tem- 
perature are taken into account. For mildly constrained T RH (e.g. T RH ^2 x 10 10 GeV), 
compatibility with the BAU observed by WMAP can be obtained for extensive new 
regions of the 62-O3 parameter space, which were previously disfavoured in the flavour- 
independent approximation. On the other hand, focusing on (normal) hierarchical light 
and heavy neutrinos, the scenario where only CP violation from the £/mns is considered 
(real R), turns out to be only marginally consistent, even for # 13 = 10°, and under mild 
bounds on T RH . Stricter T RH bounds (e.g. T RH < 5x 10 9 GeV) strongly motivate that CP 
is (also) violated in the right-handed neutrino sector. While extensive regions of the 62-63 
parameter space can produce BAU close to the WMAP range in this case, the favoured 
seesaw parameter space, clearly consistent with observations, is that of small values of 
I sin 62 1 and | sin 63 1 . 

Given the attractiveness of the mechanism of thermal leptogenesis, and the interesting 
constraints it can provide, it would be desirable to further refine the computation of the 
baryon asymmetry. Together with the expected improved bounds from LFV, electric 
dipole moments and other related observables, leptogenesis may offer valuable information 
on the right-handed neutrino masses and mixings. 
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Appendix 

A Boltzmann equations with reheating 

The efficiency factor for thermal leptogenesis introduced in Section[3]is calculated from the 
flavour-dependent Boltzmann equations [8j QUI El EE] , generalised to the MSSM case [T4"] . 
Regarding reheating, one follows the simple, yet convenient approach of Ref. [20], where 
the effects of reheating are described by a single parameter, Trh- The limitations of 
the several approaches were summarised in Section 13.1.31 In this appendix, aiming at 
completeness, we present some technical details. 

The Boltzmann equations, with z = m^/T, can be written as [T4"l [20]: 

Z— = --P+-—P+, (A.l) 

yY dY Nl 3(Z-1)X of/ /l(z) , , x (v v eqx 9 x 

ZX ~fc- = ~ ~ z kJz) ^ ' ^ 1 ~ Nl ' ' ^ ' 

ZX^ 3(Z-1)X /iW ()( eq) , , 



ZI-^ = -2e la Kz 



^ h{z) [(Y Nl - Y%) + (Y^ - Yg)\ (A.4) 



dz ^ K, 

I 3(Z-1)X K x {z) E^Va, 

+ " Y a + K a z — f 2 (z) ^ {Y Ni + Y~J. 

The above equations should be solved from z = mN 1 /T max to "infinity" (i.e. z ^> 1). 

Let us now address each of the quantities appearing in Eqs. ( 1A.1I - lA~4l) . First of all, 
let us comment on the effects of reheating. /?r and pA are the radiation energy density 
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and the energy density from the coherent oscillations of the reheating scalar field 0. The 
reheating temperature is given by 



r -li?)' (A ' 5) 

where is the decay rate of <ft, Mp\ is the Planck scale, and g* = 228.75 was already 
introduced in Eq. (TTTI) . During reheating, p^ dominates over p&. In the Boltzmann 
equations, reheating is taken into account by means of 

Z=l-^, (A.6) 
4-HPr 

which is equal to when the maximal reheat temperature T max is reached (corresponding 
to our initial conditions), and which becomes 1 after reheating. In the limit Z — > 1 (and 
X — > 1 - see definition below, in Eq. (1A.13I) ). we recover the MSSM equations without 
reheating, as given in [14] . At the maximal temperature T max , the energy density p^ 
can be calculated from the condition Z = 0, using p R = {m Nl /z) 4: 7r 2 g*/30 together with 
Eq. (IA.12I) . Notice that Eq. (IA.5[) allows to extract (appearing in Eq. (IA.1I) ) as a 



function of the reheat temperature, 



y 1 (A.7) 



Mp[ 45 

We would like to stress that in specific models of reheating after inflation, the prospects for 
leptogenesis could be significantly different. Nevertheless, this set of Boltzmann equations 
(Eqs. (1A.ll - lA~4l) ) simulates the generic constraints arising for thermal leptogenesis from 
bounds on the reheat temperature for a large class of scenarios. 

YAa are defined as Y~A a = yB/3 — YL a , with Yt a being the total (particle and sparticle) 
lepton number densities for a flavour a. Y a ^ are the densities of the (s)lepton doublets and 
Y N ft are the densities of the right-handed (s)neutrinos. The corresponding equilibrium 
number densities (in the Boltzmann approximation) are given by 

w K- eq « 1^)3 Y *\z) » Yt\ z ) w z 2 K 2 {z) , (A.8) 

with K 2 (z) (and ifi(z)) being the modified Bessel functions of the second kind. 

The matrix A, which appears in the washout term in Eq. flA.41) . is defined via the 
relation Y a = A af3 Y"a q , with Y a = Y a + Y& being the combined densities for lepton 
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and slepton doublets. Below (1 + tan 2 /3) x 10 9 GeV, where the Boltzmann equations are 
solved for the individual asymmetries Y& e , Y\ and Y& T , A is given by 



-93/110 6/55 6/55 
a mssm = | 3 / 40 _ 19 /3o 1/30 j. (A.9) 

3/40 1/30 -19/30, 

The final lepton asymmetry in each flavour is governed by Trh and by three sets 
of parameters: e^ a ,K a and K. The parameters £i )Qi , £1,5, e^ a and denote the 
asymmetries for the decays of neutrino into Higgs and lepton, neutrino into Higgsino 
and slepton, sneutrino into Higgsino and lepton, and sneutrino into Higgs and slepton, 
respectively. They are defined as 

_ T Nl e a - T Nl j a ^Nj a ~ ^Nj* a 



r~ — r~- r ~ ~ — r ~ ~ 

e ^ a ~ V (T~ + £l ' s ~~ V (T~ ~ + T~ ~ ) ' (A.1U) 

with r the decay rates of (s)neutrinos with (s)leptons in the final states. In the MSSM, 
the four decay asymmetries are equal, E\^ a = e^g = ej a = ej~, and given by Eq. fTHZ]) . 
The parameters K a control the washout processes for the asymmetry in an individual 
lepton flavour a, and K controls the source of right-handed neutrinos in the thermal 
bath. In analogy to the case without reheating, they are given by 

K a ee + r f " , K = Y. K *' K a = K%^, (A.ll) 

where Hoim^) is the "fictitious" Hubble parameter (without reheating). The latter 
is computed without taking into account p^ at T — m^, and is given by H (m Nl ) « 
l.66y/g^m% /M-pi. In the presence of p^, the Hubble parameter is modified to 



H 



87T 



1 

5 5rT < A - 12) 



and in order to match the H (real) and Hq ("fictitious") Hubble parameters, one intro- 
duces the quantity X, 



X = ( PR + /> ) . (A.13) 

PR 
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We further notice that the parameters K a are related to m ljCn introduced in Eq. ffTol) . as 
mi,« = K a m*, with m* « sin 2 (/3) x 1.58 x 10 -3 eV [H] . (A.14) 

Finally, the function f\(z) accounts for the presence of AL = 1 scatterings and ^(z) 
accounts for scatterings in the washout term of the asymmetry. They are defined as 

In 

Id + 7s,al=i = IdIi , -y + Tw,al=i = 7S/2 , (A. 15) 

where 7^ is the thermally averaged total decay rate of iV\ and 7s,al=i represents the 
rates for the AL = 1 scattering. The corresponding flavour-dependent rates for washout 
processes involving the lepton flavour a are denoted by 7^ (from inverse decays involving 
leptons £ a ) and 7w,az=i- 
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